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Abstract: We briefly describe some of our recent results for the mass spectrum and matrix 
elements using 0(a) improved fermions for quenched QCD. Where possible a comparison is made 
between improved and Wilson fermions. 



1 Introduction 



Upon discretising the QCD-Lagrangian in the Wilson formulation we find that the 
gluon piece has 0(a 2 ) errors, but for the fermion part, the additional Wilson term, 
necessary to avoid the fermion 'doubling' problem gives us 0(a) errors. Thus we 
would expect that looking at any physical quantity, for example a mass ratio, we 
have 

= ro + ari + a 2 r2 + 0(a 3 ). (1) 

m H > 

Symanzik developed a systematic perturbative progamme to 0(a n ) in which a basis 
of irrelevant operators is added to the Lagrangian to completely remove 0(a n_1 ) 
effects. Restricting this to on-shell (or physical) quantities, QXJ , enables the equations 
of motion to be used to reduce the required set of operators, in both the action and 
for improved operators in matrix elements. In this talk we briefly report our progress 



*Talk presented by R. Horsley at the 31st Ahrcnshoop Symposium on the Theory of Elementary 
Particles, September 1997, Buckow, Germany. 
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on computing the hadron spectrum and matrix elements for 0(a) improved fermions, 
where in this case we expect for a physical quantity such as a mass ratio 

!^L = ro + a ^ + 0(a 3 ). (2) 

The ground-work has been laid by the Alpha collaboration who succeeded in non- 
perturbatively calculating several improvement coefficients for (3 = 6.0/g 2 > 6.0. 
Reviews and further references may be found in JJ. We shall here neither describe 
their method nor go into details of how the mass spectrum or matrix elements are 
computed on the lattice. A description of our spectrum results is given in ref. H 
while for Wilson matrix elements see, for example, Q. Up to now we have generated 
configurations on (16 3 , 24 3 ) x 32 lattices at (3 — 6.0 and on 24 3 x 48 lattices at f3 — 6.2. 
We have used the string tension, \J~K, to determine the lattice spacing a at these 
two g 2 values. A physical value of \/K = 0.427GeV is taken to set the scale. 



2 The mass spectrum 

We have computed the p (J PC = 1 ) and nucleon (N) masses together with the a 
(0 ++ ), ax (1 ++ ) and bi (l + ~) masses. The latter three mesons being p-wave states 
are difficult to measure with our symmetric sources and our results should only 
be taken as indicative of possible trends. In Fig. [1] we plot this hadron spectrum 
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Figure 1: The hadron spectrum for light hadrons at (3 = 6.0 (circles), 6.2 (squares) using Wilson 
fermions (left picture) and using 0(a) improved fermions (right picture) against the lattice spacing, 
a. Experimental numbers are shown as stars. 

both for Wilson and 0(a) improved fermions. Roughly the same amount of CPU 
time has gone into producing each of these pictures. From eqs. ([T],|2]) we expect 
that the dominant discretisation terms are such that linear extrapolations in a for 
Wilson fermions and a 2 for 0(a) improved fermions are sufficient. Clearly with 
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only two points we must limit ourselves here to qualitative observations^. First 
we note that due to the absence of the 0(a) term for the improved fermions, the 
convergence to the continuum limit is faster. On the other hand it does appear 
that, especially for heavier particles, the signal fluctuates more. It also seems that 
at least for the p and TV particles, 0(a) improved fermions results lie closer to the 
continuum result - the Wilson results have a noticeable gradient. An alternative 
way of looking at these results is given in Fig. |2j. While it is not necessary that 
quenched QCD should reproduce exactly the continuum spectrum, in many present- 
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Figure 2: a^ 1 estimates using H = p or N as a scale normalised to the string tension for Wilson 
fermions (circles) compared to 0(a) improved masses (squares). 

day applications a hadron mass (typically the p or N) is used to set the scale. 
Potential discrepancies between these scales can be revealed by normalising them 
against the string tension (which has 0(a 2 ) errors). Again this confirms the previous 
impression: 0(a) improved fermions for the p or N perform better than Wilson 
fermions. 

3 Matrix elements 

Matrix elements, such as decay constants or moments of structure functions, are more 
complicated to calculate than masses as the operators must be appropriately renor- 
malised, so that O r = Z Q (l + boam q )0, with O = O + J2i c i a Oi. For present-day 
(3 values, non-perturbative (NP) determinations of Zq, bo and {cj} are preferable. 
(Here we are mainly interested in matrix elements in the chiral limit, so we do not 
need bo', the other irrelevant operators 0\ are only required in the improved case 
to ensure complete 0(a) cancellation.) In Fig. ^| we show Zy (V^ = q^^q) and Z A 
(A^ = ^7^755) for both Wilson and 0(a) improved fermions. We see that in all 

^Presently, for 0(a) improved fermions, most groups have concentrated on /3 = 6.0, 6.2, but 
there has been a recent interesting attempt to go to lower (3 values, 0]. 
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cases, while first order perturbation theory lies > 10% away from NP computations, 
tadpole improvement ( 77) always gives results closer to the NP line. For the Wilson 
case we expect to find larger discrepancies (of 0(a)) between the various NP esti- 
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Figure 3: Various determinations of Zy (upper row) and Za (lower row) for Wilson (left 
column) and O(a) improved fermions (right column). The tadpole improved (77) results have been 
obtained from first order perturbation results, Zo(g 2 ) = 1 + cog 2 + 0(g 4 ), 0, using the procedure 
given in 0. Also plotted are results using the local vector current between nucleon, || (filled 
squares) and p states, |i| (filled circles), and the local axial current between quark states, |Ti| 
(filled squares). Where appropriate, a simple Pade interpolation/extrapolation has been applied, 
Zo(g 2 ) = (1 + Pog 2 + log )/(l + r og 2 ) with po — r o = c o- F° r the 0(a) improved fermions we 
also give the Alpha results (filled circles) and Pade extrapolation, and from JllJ (open circle). 

mates of Z than for improved fermions (of 0(a 2 )). This seems most obvious for Zy. 
We next note that both the Alpha, ||10|| , and our non-perturbative determination of 
Zy for 0(a) improved fermions are in very good agreement with each other. From 
current conservation we expect (N\V^\N) = p^/E N x q with Xu = 2, Xd = 1- This 
at (eg) zero momentum allows a non-perturbative determination of Zy (and by), for 
both Wilson and 0(a) improved fermions. This result is exemplified in Fig. |] where 
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Figure 4: (JV|Vf jiV), left picture, (iV|Vf |AT) x E N /p x , right picture, at p = (pi,0,0) with Pl 
being the minimum possible momenta, namely pi = 27r/16a « 765MeV or 27r/24a ~ 694MeV at 
/3 = 6.0, 6.2 respectively. The expected values are Xu — 2, Xd = 1- Filled circles, squares denote 
0(a) improved fermions at /3 = 6.2, 6.0 respectively, while open squares are for fj = 6.0 Wilson 
fermions. To guide the eye, constant fits are made to the chiral limit for the right hand picture. 

we perform a check of momentum effects and restoration of rotational invariance, by 
computing (JV|Vf|iV), (JV|Vf |iV) x E N /p x . While all results for (N\V£\N) are con- 
sistent with each other, (N\Vi\N) x E^/pi shows some spread; indeed at (3 — 6.0, 
Wilson fermions seem to perform better than the 0(a) improved fermions. 

In the left picture of Fig. [| we apply the Alpha result for Za to a determination 
of gA- There seems to be an approach to the experimental value of $a (but with 
rather large 0(a 2 ) corrections). For the moments of the quark distributions in the 
nucleon, we obtain the results shown in the right hand picture in Fig. ^. The lowest 
moment, (x), derived from the operator q^^D v q is interesting, not only because it 
contains one derivative, but also because there are two distinct lattice representations 
involving diagonal or off-diagonal elements, the latter requiring for its evaluation a 
non-zero spatial momentum. (Again we always used the minimum possible.) This 
potentially allows a better study of the approach to the continuum limit. At present 
a non-perturbative renormalisation of (x) is not known; however no tadpole term 
contributes to first order perturbation theory, so that this automatically lies close 
to the TI result, so we might hope that any error in the renormalisation constant is 
small. A further problem for the 0(a) improved operators is the addition of at least 
one irrelevant operator. However, as discussed in ||, we hope that its effect is also 
small. For (x), from Fig. [|, all results appear to be reasonably consistent with each 
other, and for 0(a) improved fermions there does not seem to be any large 0(a 2 ) 
effect. However, there seems to be a larger deterioration of the signal (in comparison 
to the Wilson case) when introducing a momentum into the operator. 

In conclusion, it would seem that 0(a) Symanzik improvement does indeed bring 
us closer to the goal of calculating continuum masses and matrix elements in QCD. 
The numerical calculations were performed on the Quadrics facility at DESY-IfH. 
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Figure 5: gA — Au - Ad, left picture, and - (x n )^ moments at a scale of [i ~ 1.95GeV, 

right picture. The phenomenological results are denoted with stars (gA — 1-26 and 
by the MRS parametrisation, 0(a) improved fermion results are given by filled symbols. For 

(a;), results using a zero-momentum diagonal matrix element are given by a square while the non- 
zero momentum off-diagonal results are denoted with a circle. 
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